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The HST colours of high-redshift population III galaxies 
with strong Lya emission 



> 

in 
in 

o 



X 



Received in original form 



y—i ; Erik Zackrisson 1 *, Akio K. Inoue 2 , Claes-Erik Rydberg 1 and Florent Duval 1 

f"*"^ , Department of Astronomy, Stockholm University, Oscar Klein Center, AlbaNova, Stockholm SE-106 91, Sweden 
i 2 College of General Education, Osaka Sangyo University, 3-1-1, Nakagaito, Daito, Osaka 574-8530, Japan 

<D ■ 

' Accepted 

O 

u 
6 



ABSTRACT 

Population III galaxies, made partly or exclusively of metal-free stars, are predicted 
to exist at high redshifts and may produce very strong Lya emission. A substantial 
fraction of these Lya photons are likely absorbed in the intergalactic medium at z > 6, 
but recent simulations suggest that significant Lya emission may be detectable up 
to z w 8.5, i.e. well into the reionization epoch. Here, we argue that high-redshift 
population III galaxies with strong Lya emission can be identified in Hubble Space 
Telescope imaging data because of their unusual colours. We quantify this effect in 
some of the filters used in F-band dropout searches for galaxies at z ~ 8 and find that 
population III galaxies with high Lya fluxes may exhibit much bluer J — H colours 
at z ~ 8-10 than any normal type of galaxy at these redshifts. This colour signature 
can arise even if pop III stars account for as little as ~ 10~ 3 -10~ 2 of the stellar mass 
in these galaxies. Some of the anomalously blue objects reported in current K-band 
dropout samples do in fact meet the colour criteria for Lya-emitting population III 
galaxies. 

Key words: Galaxies: high-redshift - galaxies: stellar content - galaxies: ISM - 
galaxies: photometry 



1 INTRODUCTION 

The first population III (hereafter pop III) stars are ex- 
pected to form in isolation or in small numbers within 
~10 5 -10 6 Mp) dark matter halos at redshifts z w 10- 
60 (e.g. iTrenti fc Stiavellil |2009l ). but such stars are 
likely too faint to be detectabl e even wit h the James 
Webb Space Telescope^ (e.g. iGreif et alj 120091 ). Pop 
III stars may, however, also continue to form in the 
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it z ^ 15 (e.g. Scannapicco. Schneider & Fcrrara 


2003 


: Schneider et alj|2006l; 


Tornatore, Fcrrara & Schneider 


2007 


: 1 Johnson, Greif & Bromml 120081: 1 Johnson et al.1 2009: 


Stiavelli & Trenti 2009; Johnson 2010), and this could al- 



low their integrated light to be detected, possibly even 
with existing telescopes. Such pop III galaxieqj may 
display a number of spectral signatu res that set them 
apart from more mundane objects (e.g . Tumlinson fe Shull 
2000l:lTumlinson. Giroux^ fc Shull 2000; Schacrcr 2002. 20031 ; 
Inoue 1 2010; Raitcr, Schacrcr & Fosburv 2010; Inoue ]201ll : 
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1 http://www.jwst.nasa.gov/ 

2 here defined as galaxies which, at least over certain wavelength 
intervals, are dominated by the light from pop III stars 



IZackrisson et al1l201lh ■ In particular, the high effective tem- 
peratures (~ 10 5 K) of pop III stars imply very high fluxes 
of hydrogen-ionizing (Lyman continuum) photons, many of 
which may be converted into Lya photons at a rest wave- 
length of 1216 A after recombination in the surrounding gas. 
Under idealized conditions, this implies that pop III galaxies 
should display very high rest frame L ya equivalent wid ths 
(up to EW(Lya) « 600 - 4000 A; e.g. (Raiter et alj|20ich . 

The strength of the Lya line can admittedly be re- 
duced by a number of mechanisms. Feedback from pop III 
stars may cause Lyman continuum photons to escape di- 
rectly into the intergalactic medium (IGM), although the 
importance of this effect depends both on the pop III stel- 
lar initi al mass function (IM F) and the star formation effi- 
ciency (| Johnson et alj|2009h . Since Lya is a resonant line, 
Lya photons can also scatter repeatedly within the interstel- 
lar medium and may eventually be destroyed by dust. Cur- 
rent observations suggest that the fraction of Lya photons 
that survive increases with redshift, growing from < 0.01 
for metal-rich gala xies in the local U niverse to « 0.3 for 
galaxies at z ~ 6 llHaves et al.ll20ill ) - possibly due to a 
lower dust content in high-redshift galaxies. Since pop III 
galaxies are expected to have very little dus10, the fraction 

3 During early stages of evolution, pure pop III galaxies are ex- 
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of Lya photons that escape from such galaxies may be even 
higher. For galaxies in the reionization epoch (2 > 6), the 
Lya flux may be further reduced by absorption in the neu- 
tral 1GM. H owever, the claimed discovery of a Lya emitter 
at z Ri 8.6 jLehnert et alj|2010l ) suggests that Lya detec- 
tions are possible well into the reionization epoch, and ra- 
diative transfer simulations predict that outflows and patchy 
reionization may allow a significant fraction (w 0.1-0.5) of 
the Lya photons to evade absorpt ion up to this redshift 
(Diikstra, Mcsi nger fc Wvithel201 J ). Other teams have also 
reported a few detections of objects with high EW (Lyq) at 
z > 6 fe.g. lStark et alj|20ld : Ikashikawa et al.ll201lh . Hence, 
pop III galaxies may plausibly display strong Lya emission 
even at z > 6. This would have a pronounced effect on broad- 
band fluxes in filters tr ansmitting the Lya line, with curious 
colou rs as a result fe.g. lPello fc Schaererl2003l : lRicIiard et al.l 
2006). Here, we quantify this effect in some of the Hubble 
Space Telescope (HST) filters used in the study of z ~ 8 
galaxies. 



2 MODEL RESULTS 

To explore the impact of the Lya line on t he HST fluxes, we 
use t he Yggdrasil spectral synthesis code (| Zackrisson et al.l 
1201 ll ) w ith Starburst99 Pad o va-AGB stellar population 
spectra (jLeitherer et al ] Il999l ; IVazquez fc Leithererl 120051 ) 
for metal-enr i ched stars (pop I /II) and Schaerer ( 20021 ) and 
iRaiter et al.l l|2010h stellar population spectra for pop III 
stars. The nebular contribution to the overall spectral en- 
ergy distribution (SEP) is computed using the photoioniza- 
tion code Cloudy l|Ferland et al.lll998l ) , assuming a spherical 
geometry for the photoionized gas. The model results are 
publicly available from the lead author's homepag fl We 
have verified that the res ults are consistent with those pro- 
duced by the llnouei (|201lh model, which is based on a similar 
machinery but assumes the nebula to be plane-parallel. 

The strength of the Lya line is here regulated through 
/Lya, the fraction of Lya photons that evade both extinc- 
tion by dust within the galaxy and absorption in the neutral 
intergalactic medium (IGMjj. When computing the colours 
of pop III galaxies, we limit the discussion to /Lya values 
in the range 0-0.5, since current simulations suggest that 
/Lya > 0.5 must be e xceedingly rare for galaxies in the 
reionization epoch fe.g. iDiikstra et al.ll201ll ). The /Lya = 
case corresponds to the situation where all Lya photons are 
absorbed and there is no Lya contribution to the observed 
colours. Since we are here focusing on galaxies in the reion- 
ization epoch, the flux at wavelengths shortward of Lya is 
always assumed to be zero due to IGM absorption. 

Concerning the possibi lity of Lyman contin uum leakage 
directly into the IGM (e.g. IJohnson et al1l2009l ). we will for 
simplicity consider only the limiting cases: either complete 
leakage of ionizing photons (/ CS c,LyC = 1, implying a purely 
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Figure 1. The temporal evolution of the rest frame Lya equiv- 
alent widths for instantaneous-burst populations with various 
metallicities and IMFs (see main text for details). To allow a com- 
parison to previous studies, complete transmission of Lya photons 
is assumed (i.e. /Lya = !)■ The abrupt end of the cyan line after 
t ?3 3 Myr signals the death of the least massive stars (50 Mg) 
in the pop III. 1 scenario. During the first few Myr, the pop III 
models are able to reach EW(Lya) Ri 1000-3000 A, whereas pop 
I and II models are confined to EW(Lya) ^ 500 A. 



stellar SED) or no leakage whatsoever (/ eS c,LyC = 0, which 
maximizes the nebular contribution to the SED). 



2.1 Lya equivalent widths 

In Fig. [T] we plot the predicted evolution of the intrinsic Lya 
rest-frame equivalent width EW(Lya) for various metallic- 
ities and IMFs, under the assumption of an instantaneous 
burst of star formation and fh ya = 1. As in I Zackrisson et all 
l|201ll) . the pop I (Z = 0.02 0) and pop II {Z = 0.0004) mod- 
els assume a lKroupal ()200ll ) IMF, representative for star for- 
mation in the local Universe. In the case of pop III, three 
different IMFs are considered: an extremely top-heavy IMF 
(pop HI.l) with power-law slope a = 2.35 (AN /AM oc M~ a ) 
throughout the mass range 50-500 Mq; a moderately top- 
heavy, log-normal IMF (pop III. 2) with characteristic mass 
M c = IOMq , dispersion a = 1 Mq and tai l s exte nding 
from 1-500 Mq; and finally the same iKroupal <|200ll ) IMF 
as adopted for pop I /II. 

At ages up to a few Myr, the EW(Lya) of the pop III 
models are in the range w 1000-3000 A, which is sufficient 
to have a pronounced effect on HST broadband fluxes of pop 
III galaxies, even if just a minor fraction (/Lya > 0.1) of the 
Lya photons are transmitted to the observer. By contrast, 
the pop II model peaks at EW(Lya) ~ 500 A and the pop I 
model at ~ 200 A. These EW(Ly a) predictions are la rgely 
consistent with those presented bv lRaiter et al.l l|2010l) . 



pected to have no dust, but Lya photons originating from hybrid 
galaxies consisting of both pop III and pop II/I stars may well be 
subject to extinction effects 

4 Yggdrasil model results available at: www.astro.su.se/~ez 

5 Our /Lya should not be confused with the usual definition of 
the Lya escape fraction, as /Lya is the product of Lya escape 
fraction and the IGM transmission 



2.2 The colour signatures of pop III galaxies with 
strong Lya emission 

In Fig. [2] we plot the predicted evolution of the HST/WFC3 
Jiih — Hieo colour for pop I, II and III galaxies at z = 8.5 
with and without residual Lya emission. In the absence of 
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Figure 2. Predicted HST/WFC3 J125 — J^i60 colours at z = 8.5. a) The predicted evolution for pop III galaxies with 
IMFs and different degrees of nebular and Lya emission, b) The predicted evolution for pop I, II and III galaxies with 
IMFs and nebular emission but varying of Lya fractions. The bluest models at this redshift are the pop III galaxies with high Lya 
fractions (/Lya = 0.5; blue solid line in both panels). 



any significant Lya contribution in the J125 filter, nebular 
emission (red line in Fig. [5^.) makes the J125 — -Hieo colour 
redder compared to a purely stellar SED (orange line in 
Fig. However, once a non-zero fhya is assumed (green, 
cyan and blue lines for fhya = 0.1, 0.3 and 0.5), J125 — 
Hieo become progressively bluer. A similar result was also 
presented bv lPello fc Schaererl ([20030 - 

As shown in Fig.[2p, young pop III galaxies with strong 
Lya emission are moreover predicted to be bluer than both 
pop I and II galaxies with similar fhya at this redshift. As 
dust extinction likely renders fhya lower in pop I /II galaxies 
and further reddens the colours, it therefore seems reason- 
able to search for population III galaxies among the objects 
with the very bluest J125 — -H16O colours. 

The pop III models in Fig. [2] are all based on a lKroupal 

([20011 ) IMF. Since more top-heavy IMFs can produce higher 
EW(Lya) (see Fig. [TJ, one would naively expect that such 
pop III models should give rise to even more extreme colours. 
However, we find this not to be the case. While the pop III 
IMF is important for determining the duration of the lumi- 
nous ph ase of a pop HI galaxy, the colours of pop III galax- 
ies with iKroupa) (|2001j ). pop III.2 and pop III.l IMFs are 
nearly identical at young ages, despite different EW(Lya). 
This is because the higher Lyman continuum flux associated 
with a more top-heavy IMF boosts the strength of the Lya 
emission line and the relative strength of nebular continuum 
directly longward of Lya. In these filters, the two effects bal- 
ance each other to make the overall colours fairly insensitive 
to the overall IMF. 

In Fig. we show the predicted redshift evolution of 
the J125 — -H16O colour for newborn (1 Myr old) galaxies 
of various metallicities and fhya (the same models as in 
Fig. [2p). Throughout the redshift interval 2 ~ 8-10, pop 
III galaxies with significant Lya emission (/Lya = 0.5) dis- 
play J125— -ffi60 ^ —0.6. Under the assumption that no other 
galaxies display /Lya > 0.5, pop III objects are the only nor- 
mal reionization-epoch (2 > 6) galaxies that are expected 



to display such colours. Objects with J125 — Hmo —0.6 
(gray region in Fig. therefore represent excellent popula- 
tion III galaxy candidates. An age of 1 Myr is adopted in 
this plot since this produces the bluest colours possible for 
pop III, II and I galaxies (see Fig. [2p). Higher ages would 
imply redder colours, as demonstrated in Fig. [2] for an in- 
stantaneous burst, but pop III galaxies with constant star 
formation rates can in principle retain their very blue colours 
( J125 - #160 < -0.6) for up to ~ 10 s yr. 

Intriguingly, reionization-epoch objects with colours ex- 
pected for young pop III galaxi es with Lya emission may al- 
ready have been detected. The lTaniguchi. Shiova fc Trumpl 
(2010) compilation of 2 « 8 objects contains one objects at 
J125 — ffieo < —0.6 and one at J125 — #160 < —0.5. In Fig.0] 
we compare the predicted J125 — #160 and Y105 — J110 colours 
of Lyq-emitting pop H I galaxies at 2 ~ 8 to those of the 
iTaniguchi et al.l (|2010h sample and the Yios-band dropout 
criteria (light gray region) typic ally used to select z ~ 8 tar- 
gets from multiband HST data |f3ouwens et al.|[20Tll ). This 
diagram also includes a dark gray, pop III galaxy zone. The 
extent of this area has been derived by studying the evolu- 
tion of all pop I, II, III galaxies at z = 7-9 in our model 
grid and identifying the region in which only pop III galax- 
ies with strong Lya emission should appear. Indeed, the two 
bluest objects (red markers with black rings) in J125 — i/160 
are broadly consistent with the colour criteria for pop III 
galaxies, at a likely redshift of 2 w 8.1-8.9. While not in- 
cluded in this plot, we note that the sample of iMcLure et al.l 
i|201ll) includes two more objects with J125 — Hkm < —0.5 
(HUDF-2641 and ERS-9041, in their notation). 

iDunlop et all (|201 ll ) and IWilkins etahl (|201ll ) caution 
that spurious detections of anomalously blue objects such as 
these are expected close to the detection thresholds. Follow- 
up spectroscopy is therefore required to unravel the exact 
nature of these pop III candidates. If their colours are ap- 
proximately correct and indeed due to Lya, our models pre- 
dict these objects to display rest-frame EW(Lya) > 200 A. 
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Figure 3. The redshift evolution of the HST/WFC3 J125 - H 1(i0 
colour for 1 Myr old galaxies with various degrees of nebular 
and Lya emission. The models used are the same as in Fig. l2b. 
Among reionization-epoch galaxies, only pop III galaxies with 
strong Lya emission (blue solid line) display J125 — #160 ^ —0.6, 
which makes objects in this colour range (gray region) excellent 
pop III candidates throughout the redshift interval z fa 8—10. 

Fig.[4]also reveals that pop III galaxies with strong Lya 
emission require higher redshifts than normal galaxie s to en - 
ter a given dropout-selected sample. iBouwens et al.l (|201ll ) 
adopt a Yios-band dropout criteria of Y105 — J125 > 0.45 and 
J105 — Hieo < 0.5 for their z ~ 8 sample (light gray region), 
and estimate that this will select objects in the z « 7-9.2 
interval. However, since the Yios and J125 filters partly over- 
lap, a strong Lya line at z « 8 can affect both fluxes and 
prevent the Y105 — J125 > 0.45 condition from applying. Our 
models for young pop III, /Lya = 0.5 galaxies suggest that 
such objects can escape the dropout criterion until z > 8.1 
and then remain a Vios-band dropout (without becoming a 
J125 band dropout) as long as z ^ 1 0.4. At 7.0 < z < 8.1 , 
these objects would instead enter the IBouwens et alJ (|201ll ) 
Z850-band dropout sample (assumed to contain objects in 
the range z « 6-7.5). Hence, dropout selection schemes of 
this type may lead the redshifts of strong Lya-emitters to 
be underestimated by up to A(z) « 0.5-1. 



3 DISCUSSION 

The results in Sect. 12.21 are valid for pure pop III galax- 
ies, i.e. objects consisting of pop III stars only. However, 
the detection of such galaxies at z ~ 8.5 in the deepest 
HST fields (J125 ^ 29.9) would require stellar population 
masses of > 10 5 M Q for a pop III.l IMF and > 10 6 Af Q for a 
Kroupal [2 00 1) IM F, ass uming fhya ^ 0-5. In the models of 
Trenti et all (|2009l ) and IStiavelli fc Trentil (|2009l ). the most 
massive pure pop III galaxies at z > 7 have halo masses of 
M ~ 10 8 Mq and baryonic masses of M ~ 10 7 M Q . Hence, 
the star formation efficiency (the fraction of the total halo 
gas mass converted into stars) would need to be very high 
(e > 0.01-0.1) to make such objects detectable in current 
data. 




Figure 4. The location of z ~ 8 pop III galaxies in the 
J125 — H16O vs - Yios ~ J125 plane. Blue lines represent pop III, 
instantaneous burst models with /Lya = 0.5 at 2 = 7.5, 8.0, 8.5 
and 8.9. Along these tracks, triangles, circles and squares indi- 
cate ages of 1, 10 and 100 Myr. Th e light gray reg i on in dicates 
the Y"-band dropout criteria used by|Bouwcns et al. (2011) to se- 
lect z ft: 8 galaxies. The dark gray region represents the region 
expected to be populated exclusively by pop III galaxies with 
strong Lya emission at z = 7—9. Red markers represent the ob- 
jects in the Taniguchi et al. compilation of z ft 8 galaxies. 
Two objects in this sample (red symbols surrounded by black cir- 
cles) are potential po p II I gal axy candida tes. The red star marks 
the position of the Leh nert et alJ l|201Cjl spectroscopically con- 
firmed z ft: 8.6 Lya-emitter. To indicate th e typical colours of ob - 
jects in deep HST fields, the objects in the Cam eron et alJ l|2010T ) 
HUDF09 source catalog have been included as black dots. 



Galaxies in which pop III and pop I/II stars form in 
parallel can attain much higher total masses, but any pop 
III spectral signature may at the same time be d iluted be- 
yond recognition. While lJimenez &: Haima n (2006) interpret 
the Lya properties of galaxies at z » 3-4 as due to a very 
large mass fraction (> 0.1) in pop III stars, simulations tend 
to fav our a much smaller contribution (e.g. ISalvaterra et al.l 
2011). To explore the requirements for a hybrid galaxy to 
display the pop III spectral signatures discussed in this pa- 
per, we define a mass ratio / pop m of pop III stars within 
a galaxy: / pop m = M pop m/(M pop m + A/ pop I/n ), where 
the different masses represent the gas masses converted into 
stars prior to the age at which the galaxy is observed. Since 
a young burst of pop III stars has a rest-frame UV stellar 
M/L ratio that is much smaller than that of its pop I/II 
host galaxy, the pop III component can dominate the UV 
flux even if its mass is significantly lower. 

Obviously, the / pop in ratio required to produce ex- 
tremely blue colours due to Lya domination depends on 
many factors, including the pop III IMF and the star for- 
mation history of the pop I/II component. However, assum- 
ing a 1 Myr old pop III burst located in a 300 Myr old pop 
I galaxy with a constant star formation rate (both having 
/Lya = 0.5) at z — 8.5, we find that extremely blue colours 
{J105 — Hieo ^ —0.6) can be produced fo r / pop m ^ 7x l0~ 2 
if the pop III IMF is similar to that of iKroupal j|200ll ), for 
/ pop in ^ 1 x 10 -2 in the case of a pop III.2 IMF and for 
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/pop in ^ 4 x 1CP 3 in the case of a pop III. 1 IMF. If the pop 
III burst takes place in a more passively evolving pop I /II 
galaxy (modelled as in instantaneous burst ), the ra t io can 
be as low as / pop m ~ 5 x 10 -3 even for a iKroupal (|200lh 
pop III IMF and « 3 x 10~ 4 for a pop III.l IMF. 

The models presented in this paper assume the stel- 
lar IMF to be fully sampled. This assumption breaks 
down in systems containing small numbers of stars, with 
drastically different colour p redictions as a result (e.g. 
ICervino fc Valls-Gabaudlinol ). While pure pop III galaxies 
are expected to be low-mass systems with limited numbers 
of stars, IMF sampling effects are unlikely to jeoparize colour 
signatures dominated by strong Lya fluxes. As discussed in 
Sect. 12.21 our colour predictions are very similar for the three 
pop III IMFs considered, and this holds as long as nebular 
emission dominates over direct star light. To have any signifi- 
cant impact on the colours, IMF sampling effects would have 
to result in a situation where almost no pop III stars with 
mass > 10 Mq are present, as this would imply very little 
ionizing radiation. This is impossible for the two top-heavy 
IMFs consid ered (pop III.l and III. 2), but may in principle 
occur for the lKroupal q2001r > IMF. However, this would also 
dramatically reduce the luminosity of such systems. We es- 
timate that a stellar population mass of > 10 7 " 8 Mq would 
be required to make such pop III galaxies detectable in cur- 
rent HST data at«» 8.5. This is orders of magnitude higher 
than the mass scale where I MF sampling is likely to affect 
photoionization calculations (|Villaverde et al.ll2010l ). Hence, 
for pop III galaxies in the luminosity range relevant for cur- 
rent observations, IMF sampling is not likely to be an issue. 

Finally, we note that while no normal pop I/II galax- 
ies are expected to to enter the Jno — -ffi25 < —0.6 pop 
III region at z > 5 (see Fig. [3|, it is conceivable that 
other unusual high-redshi ft objects cou ld display similar 
colours. Gas cooling (e.g. iDiikstral 120091 ) or fast accretion 
shocks l|Dopita et alj|2011f ) in massive halos may also pro- 
duce strong Lya emission with high equivalent widths, and 
so can accreting b lack holes (e.g. lHaiman fe Reea l200ll: 
Ijohnson et al.l l2011). Since the identification of any of these 
mechanisms at z > 8 would be a remarkable discovery in 
its own right, follow-up studies of objects with pop III- like 
colour signatures are definitely justified. 
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